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Abstract. 

Powerlaws have long been used to describe the spectral dependence of aerosol 
extinction, and the wavelength exponent of the aerosol extinction powerlaw 
is commonly referred to as the Angstrom exponent. The Angstrom exponent 
is often used as a qualitative indicator of aerosol particle size, with values 
greater than two indicating small particles associated with combustion byprod- 
ucts, and values less than one indicating large particles like sea salt and dust. 

In this study, we investigate the relationship between the Angstrom expo- 
nent and the mode parameters of bimodal aerosol size distributions using 
Mie theory calculations and Aerosol Robotic Network (AERONET) retrievals. 

We End that Angstrom exponents based upon seven wavelengths (0.34, 0.38, 

0.44, 0.5, 0.67, 0.87, and 1.02 /irn) are sensitive to the volume fraction of aerosols 
with radii less then 0.6 /irn, but not to the hne mode effective radius. The 
Angstrom exponent is also known to vary with wavelength, which is com- 
monly referred to as curvature; we show how the spectral curvature can pro- 
vide additional information about aerosol size distributions for intermedi- 
ate values of the Angstrom exponent. Curvature also has a significant effect 
on the conclusions that can be drawn about two-wavelength Angstrom ex- 
ponents; long wavelengths (0.67, 0.87 /irn) are sensitive to hne mode volume 
fraction of aerosols but not hne mode effective radius, while short wavelengths 
(0.38, 0.44 /irn) are sensitive to the hne mode effective radius but not the 
hne mode volume fraction. 
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1. Introduction 

Knowledge of the aerosol optical thickness (AOT) throughout much of the shortwave 
spectral region (~0.3-5 /irn) is necessary to compute the shortwave aerosol radiative forc- 
ing at the surface and the top of the atmosphere. The AOT is easily measured in discrete 
spectral intervals with sunphotometers located at the surface, but gas and water vapor 
absorption prevent the measurement of AOT at all wavelengths of interest. This diffi- 
culty is easily circumvented because Angstrom [1929] noted that the spectral dependence 
of extinction by particles may be approximated as a powerlaw relationship: 

r(A) = TiA" a , (1) 

where r(A) is the aerosol optical thickness (AOT) at the wavelength A, T\ is the approx- 
imated AOT at a wavelength of 1 /mi sometimes called the turbidity [ Angstrom , 1964], 
and a has come to be widely known as the Angstrom exponent. 

In addition to being a useful tool for extrapolating AOT throughout the shortwave 
spectral region, the value of the Angstrom exponent is also a qualitative indicator of 
aerosol particle size [ Angstrom , 1929]; values of a < 1 indicate size distributions dominated 
by coarse mode aerosols (radii > 0.5 /mi) that are typically associated with dust and sea 
salt, and a > 2 indicating size distributions dominated by foie mode aerosols (radii 
< 0.5 /irn) that are usually associated with urban pollution and biomass burning [Eck 
et al, 1999; Westphal and Toon, 1991]. Kaufman et al. [1994] demonstrated that the 
Angstrom exponent can be a good indicator of the fraction of small particles with radii 
r = 0.057-0.21 /mi relative to larger particles with radii r = 1.8-4 /rm for tropospheric 
aerosols. 
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Since the Angstrom exponent is easily measured using automated surface sunphotom- 
etry [ Holben et al, 1998] and is becoming increasingly accessible to satellite retrievals 
[Nakajima and Higurashi, 1998; Higurashi and Nakajima, 1999; Deuze et al., 2000; Ig- 
natov and Stowe, 2002; Jeong et al, 2005], the true utility of this parameter lies in its 
empirical relationship to the aerosol size distribution. For instance, a has been used to 
characterize the maritime aerosol component at island sites [ Kaufman et al, 2001; Smirnov 
et al, 2002, 2003], biomass burning aerosols in South America and Africa [Dubovik et al, 
1998; Reid et al, 1999; Eck. et al, 2001b, 2003], and urban aerosols in Asia [Eck. et al, 
2001a]. Measurements also indicate that the Angstrom exponent varies with wavelength, 
and that the spectral curvature of the Angstrom exponent contains useful information 
about the aerosol size distribution [King and Byrne, 1976; King et al, 1978; Eck et al, 
1999; Eck. et al, 2001a, b, 2003; Kaufman, 1993; O'Neill and Royer, 1993; O'Neill et al, 
2001a, b, 2003; Villevalde et al, 1994], 

The focus of this paper is to explore the relationship between the spectral dependence 
of extinction and the size distribution of atmospheric aerosols. We begin by illustrating 
the sensitivity of a to the median radius of monomodal aerosol size distributions, using 
multi- wavelength Mie calculations of r(A) for 38 monomodal lognormal aerosol size dis- 
tributions. Then we explore the relationship between a and 45 bimodal lognormal aerosol 
size distributions, demonstrating that a is more sensitive to the fine mode volume fraction 
than the fine mode median radius. Next, we apply the same technique to explore the in- 
formation content in the wavelength-dependence of a (i.e., the curvature of a). Finally, we 
discuss application of the Angstrom exponent and the spectral curvature for setting limits 
on the possible size distributions associated with aerosol optical depth measurements. 
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Before proceeding further, however, a word about the origination of the term “Angstrom 
exponent” is in order. The term “Angstrom exponent” originates from an early treatise 
by Anders Angstrom that provides the first documentation of Equation (1) currently 
available in english [ Angstrom , 1929]. However, that article cites an even earlier laboratory 
study whereby Lundholm documented the powerlaw relationship for the absorption of 
thin powders at Uppsala in 1912, prompting at least one author to point out that we 
are honoring the wrong person [ Bohren , 1989]. Unfortunately, the Lundholm citation in 
Angstrom [1929] is obscure and perhaps incorrect, as our library staff was unable to locate 
it. Compounding matters even further, it is not even clear whether Lundholm or Lindholm 
documented the powerlaw in 1912, as Angstrom [1929] uses both forms of spelling. Our 
library staff did find references to a F. Lindholm at Uppsala in Science Abstracts, but none 
of the abstracts mention the spectral dependence of particulate extinction. Subsequent 
references to Lundholm and Lindholm on this topic have not appeared in the atmospheric 
literature, and Angstrom claimed full credit for introducing the methods of evaluating 
atmospheric turbidity parameters in later articles [ Angstrom , 1930, 1961, 1964], So who 
should get the credit for originating and promoting Equation (1)? After presenting the 
discussion above, we maintain the standard nomenclature. Angstrom admittedly did not 
originate the empirical powerlaw, but he did publish at least four articles documenting the 
relationship between a and particle size; meanwhile, we are unable to obtain any public 
documentation of Equation (1) by Lundholm or Lindholm, making it impossible for us to 
fairly evaluate this person’s contribution. 
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2. Extinction Calculations 

Atmospheric aerosols are never limited to a single particle size, so we seek realistic 
polydisperse aerosol size distributions for our calculated values of r. Although Angstrom 
[1929] documented the relationship between the Angstrom exponent and single particle 
sizes, Junge [1955] was the first to explore a relationship between the Angstrom exponent 
and polydisperse aerosol size distributions. By assuming a powerlaw relationship for the 
number density (. N ) of aerosols as a function of radius (r): 


dN(r) 


= Ar~ 


( 2 ) 


d log(r 

where A and v are the coefficients that characterize the size distribution, Junge [1955] was 
able to show that u « a + 2 for nonabsorbing aerosols with a > 1 (see also Junge [1963]). 
The simplicity of this relationship has found widespread appeal, and the “Junge distri- 
bution” has subsequently appeared in numerous textbooks [ McCartney , 1976; Measures, 
1984; Stephens , 1994; Seinfeld and Pandis , 1998; Pruppacher and Klett, 1997]. However, 
the simple relationship v « a + 2 does not hold for small values of a (a < 1) or for 
absorbing aerosols [Junge, 1955; Tomasi et al, 1983; Cachorro et al, 1993; Cachorro and 
de Frutos, 1995], and the Junge powerlaw distribution is mainly of historical interest. 

Subsequent research has shown that Equation (2) does not accurately describe atmo- 
spheric aerosol size distributions, and that multimodal lognormal distributions are more 
appropriate [Davies, 1974; Whitby, 1978; Ott, 1990]: 


dV{r) ^ Ci 

- 1 L ex P 


(lnr — lni?j ) 5 

2of 


( 3 ) 


d In r ’ 

where Ci represents the particle volume concentration, Ri is the median or geometric 
mean radius, a* is the variance or width of each mode, and n is the number of lognormal 
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aerosol modes. Here, we have switched to the volume concentration dV{r ) = 4/37T r 3 d,N(r) 
because the optical effects of atmospheric aerosols are more closely related to their volume 
than their number [Whitby, 1978; Seinfeld and Pandis, 1998]. Whitby [1978] observed 
thousands of aerosol size distributions, and found that most could be described with 
three volumetric modes: a nuclei mode with geometric mean radii of 0.0075-0.020 /m i, 
an accumulation or fine mode with geometric mean radii of 0.075-0.25 //rn , and a coarse 
particle mode with geometric mean radii of 2.5-15 fii n. An additional mode with a 
geometric mean radius of ~0.5 /mi was observed by Kaufman et al. [1994] after the Mt. 
Pinatubo volcanic eruption in June, 1991; they attributed this mode to the stratospheric 
aerosol loading at that time. 

We can use Equation (3) to calculate the extinction of a polydisperse distribution of 
aerosols with refractive index m in an atmospheric column of height Z: 


r(A) = 


3Qext('m, r, A) dV 

4 r d In r 


d In r dZ , 


where Q ex t is the extinction efficiency, which we calculate exactly for spherical aerosols 
using Mie theory [ Wiscombe , 1980]. The extinction is sensitive to the size of the particles 
in the distribution, which we illustrate at two wavelengths in Figure 1 for monomodal size 
distributions with variable median radii, but identical mass and width. (Here, we have 
plotted the results in terms of the mean effective radius ( R e ff ) of the size distribution, 
which we define below). Note that the fine mode particles (R e ff < 0.25 /mi) have a 
much greater effect on the AOT at the visible wavelength (0.5 /mi) than at the near 
infrared wavelength (1.0 /mi) . Likewise, the coarse mode particles ( R e ff > 2 /mi) provide 
similar contributions to the AOT at both wavelengths in this example. The nuclei mode 
(0.0075 /mi < R e ff < 0.020 /mi) contributes very little to the AOT at both wavelengths, 
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so we do not calculate the optical effect of this mode for the remainder of this article. We 
show in later sections that the interplay of the aerosol optical depth at the visible and 
near infrared wavelengths with the fine and coarse mode particles provides information 
about the aerosol size distribution. 

We use the effective radius throughout this article, which is defined as the average radius 
weighted with the geometrical cross-sectional area [ Hansen and Travis, 1974], 

= £m^lnr 
£ff r-^dhrr’ 

and is often used for calculations involving droplet size distributions. This is because 

spherical drops scatter radiation in proportion to their geometrical cross-sectional area 

(7rr 2 ), so the single-scatter properties of size distributions are more closely related to 

the effective radius than to the median radius [ Hansen and Travis, 1974; Liou, 1992; 

Mishchenko et al, 1997]. This is useful because the nomenclature for describing lognormal 

aerosol size distributions is not universal, and size distributions with the same effective 

radius and the same effective variance will have similar single-scatter albedos [Mishchenko 

et al, 1997]. The relationship between the effective radius and the median radius for 

monomodal size distributions (n = 1 in Equation 3) is shown for three mode widths in 

Figure 2. 

3. Spectral dependence of r 

In practical applications, the extinction of atmospheric particles are measured at two 
or more wavelengths in the ultraviolet, visible, or near infrared spectral regions and the 
Angstrom exponent is calculated from the slope of the linear regression of the logarithm 
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of Equation (1), viz: 

lnr(Aj) = lnri — a In X tl (6) 

where the subscript i emphasizes that the measurements are obtained at discrete wave- 
lengths. The measurement wavelengths are chosen to avoid spectral regions with signifi- 
cant gas absorption, although spectral regions with some ozone absorption may be used 
by including a compensating adjustment. Measurements are typically obtained with nar- 
rowband wavelength filters that are about 10 nm wide and span the 0. 3-1.1 /m i spectral 
response region of silicone photodetectors [Shaw et al, 1973]. 

The Angstrom exponent itself varies with wavelength, and a more precise empirical re- 
lationship between aerosol extinction and wavelength is obtained with a 2nd-order poly- 
nomial [King and Byrne , 1976; Eck et al, 1999; Eck. et al., 2001a, b, 2003; Kaufman, 
1993; O'Neill et al, 2001a, 2003]: 

lnr(Aj) = a 0 + an In A* + a 2 (ln A*) 2 . (7) 

Here, the coefficient <22 accounts for a “curvature” often observed in sunphotometry mea- 
surements. Note that for the special case of <22 = 0, the coefficients ao and oj equate to 
a 0 = lnri and oq = —cc. 

Some authors have noted that the curvature is also an indicator of the aerosol particle 
size, with negative curvature indicating aerosol size distributions dominated by the fine 
mode and positive curvature indicating size distributions with a significant coarse mode 
contribution [ Kaufman , 1993; Eck et al, 1999; Eck. et al, 2001b; Reid et al, 1999]. This 
is illustrated graphically in Figure 3, which shows the spectral dependence of extinction 
for two plausible aerosol size distributions that have the same Angstrom exponent of 
a — 2. The AOT with negative curvature was calculated for a monomodal size distribution 
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with a median radius of R = 0.12/un to characterize a fine mode distribution, and the 
AOT with positive curvature was calculated for a bimodal aerosol size distribution with 
a significant coarse component of C crs / (Cf ine + C crs ) = 0.4. The negative curvature 
associated with the monomodal aerosol size distribution is a natural byproduct of the 
Mie calculation, and is typical of the accumulation mode sizes common to atmospheric 
aerosols. Positive curvature occurs for the bimodal distribution because the coarse mode 
AOT at the near-infrared wavelengths is a much greater fraction of the visible AOT than 
the same ratio induced by fine mode aerosols. That is, AOT(A = 1 /nn)/AOT(A = 0.5 /im) 
at R e ff = 2 /im is much greater than AOT(A = 1 /tm)/AOT(A = 0.5 /im) at _R e // = 0.1 /im 
(see Figure 1, for example). We discuss this further in Sections 4 and 5. 

4. Response of a. to changes in polydisperse aerosol parameters 

We study the impact of particle size on the spectral variability of extinction by calculat- 
ing r(A) for a variety of plausible aerosol size distributions (via Equation 4). That is, we 
test the sensitivity of the empirical coefficients a, cq, and a 2 (in Equations 6 and 7) to the 
variability of the parameters of bimodal aerosol size distributions Ci, Ri, cti (in Equation 
3), using exact calculations of r(A) at seven wavelengths (0.34, 0.38, 0.44, 0.5, 0.67, 0.87, 
and 1.02 /tn a). We normalize the aerosol optical depth at A = 0.44 /im by scaling 
but note that our computations of a, cq, and a 2 are independent of this scale factor. This 
would not be the case for real atmospheric measurements, where the determination of 
these empirical coefficients at low optical depths could be altered because of instrument 
accuracy and resolution. 
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4.1. Monomodal aerosol size distributions 

We begin by calculating the Angstrom exponent (a) for a set of spherical monomodal 
aerosol size distributions with a variety of median radii (Ri = 0.04 to 7 /an). We chose two 
widths for our monomodal size distributions consistent with the fine and coarse modes 
in the climatology of Dubovik et al. [2002]: o = 0.38 for median radii Ri < 0.6 /tm, and 
o = 0.75 for median radii Ri > 0.6 /im. We simulated an internal aerosol mixture for our 
distributions, tuning the volume fractions of water, ammonium sulfate, and black carbon 
until we achieved a minimum y 2 fit to a refractive index of m = 1.37 — 0.003i at four 
wavelengths: 0.44, 0.67, 0.87, and 1.02 /im [ Schuster et al., 2005]. This technique results 
in a plausible aerosol mixture with a spectrally variable refractive index. We use linear 
regression of Equation (6) to calculate the associated empirical Angstrom exponents. The 
sensitivity of the Angstrom exponent to the effective radius of the size distribution is 
shown in Figure 4a. Although we show two different mode widths (a) for the small and 
large effective radii in this test in order to maintain values similar to the Dubovik et al. 
[2002] climatology, the results are qualitatively the same for either mode width at all 
effective radii; that is, a marked decrease in the Angstrom exponent for particles smaller 
than about 0.6 /tm and relatively little sensitivity to larger particles [ Cachorro et al., 
2000], 

Figure 4a indicates that the Angstrom exponent is quite sensitive to the effective radius 
of size distributions that are consistent with fine mode aerosols (R e ff < 0.25 /tm), but 
provides little information about the effective radius of coarse monomodal aerosols (as 
evidenced by the nearly constant value shown for R e jf > 2 /tm). This is because the spec- 
tral variability of extinction diminishes for particles larger than the incident wavelength 
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(i.e., the visible and near infrared AOTs in Figure 1 are nearly equivalent for coarse mode 
particles). Nonetheless, the low a associated with coarse mode aerosols can significantly 
reduce the Angstrom exponent when they are included in bimodal aerosol size distribu- 
tions. We explore the impact of coarse mode aerosols on the Angstrom exponent for 
bimodal aerosol size distributions in Section 4.2. 

We also show the spectral curvature ( 02 ) as a function of the effective radius for 
monomodal aerosol size distributions in Figure 4b. Note that the curvature is always 
negative for size distributions typical of fine mode aerosols (R e ff < 0.25 /im). The curva- 
ture can become positive at larger radii, but the transition radius from negative to positive 
curvature is dependent upon the width (a) of the distribution. We explore the impact 
of coarse mode aerosols on the curvature of bimodal aerosol size distributions further in 
Section 5. 

4.2. Bimodal aerosol size distributions 

Atmospheric aerosols are rarely monomodal [ Whitby , 1978], so Figure 4a is of limited 
value for understanding the relationship between atmospheric aerosol size distributions 
and the Angstrom exponent. Hence, we repeated our calculations, this time limiting 
the fine mode median radius from 0.06 to 0.3 /irn and including a coarse mode with 
R crs = 3.2 //rn and a crs = 0.75. We also varied the ratio of the concentration of aerosols in 
the fine mode to the combined aerosol concentration [C f ine / (C f ine + C crs )\ from 0.1 to 1, 
scaling both modes as necessary to maintain a normalized optical thickness at A = 0.44 /irn. 
The results are shown in Figure 5, where we have defined the fine mode effective radius as 
the effective radius for the distribution of particles with radii less than 0.6 /irn (the mode 
separation radius recommended by AERONET at http://aeronet.gsfc.nasa.gov). 
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The uppermost line in Figure 5 represents a monomodal aerosol size distribution and 
is a subset of the line in Figure 4a with a = 0.38. The lines below the uppermost line 
in Figure 5 represent the influence of including various concentrations of coarse mode 
aerosols. The figure indicates that increasing the concentration of coarse mode particles 
(i.e., C fi ne / C total I) reduces the Angstrom exponent and dampens its sensitivity to the 
fine mode effective radius. This latter point is evidenced by the reduction in the slope of 
the C fine/C 't ota i lines with decreasing C f ine / C tota i at fine mode effective radii common to 
atmospheric aerosols (i.e., 0.1 < R e ff(fine) < 0.25 /mi). This reduction in slope occurs 
because we are adding coarse mode particles with spectrally flat extinctions to our system 
of particles, reducing the overall spectral variability . 

ft is also evident in Figure 5 that an increase in the Angstrom exponent does not neces- 
sarily correspond to a decrease in the fine mode effective radius. For instance, a bimodal 
size distribution with C fi ne /C to tai = 0.8 and R e ff(fine) = 0.07 /ini has an Angstrom ex- 
ponent of a = 2.1 in Figure 5, but the same modal concentration ratio and a fine mode 
effective radius of R, e j f(fine) = 0.09 /mi has an Angstrom exponent of 2.4. This is con- 
trary to the conventional wisdom of Figure 4, which holds that the Angstrom exponent 
decreases with increasing particle size (all else being equal). This “unconventional” be- 
havior shown in Figure 5 occurs because the coarse mode exhibits a strong influence on 
the Angstrom exponent when added to distributions composed of the least efficient parti- 
cles (i.e., the smallest ones). The fine mode extinction becomes more efficient as the fine 
mode effective radius increases beyond the smallest values shown in Figure 5 , however, 
and the conventional Angstrom exponent behavior appears once again (i.e., a decreases 
as Ref f increases). 
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5. Response of spectral curvature to changes in polydisperse aerosol 
parameters 

Recall that the traditional powerlaw for aerosol extinction is strictly an empirical re- 
lationship (Equation 1), and that the wavelength dependence of aerosol extinction is 
more accurately described using a 2nd-order polynomial (Equation 7). Hence, we repeat 
our study of the wavelength-dependence of bimodal lognormal aerosol size distributions 
by once again calculating the extinction of spherical aerosols at seven wavelengths (via 
Equations 3 and 4), this time determining the empirical coefficients a 0 , Gq, and a 2 of 
Equation (7) through polynomial regression. We tested bimodal size distributions with 
all permutations of mode sizes, widths, concentrations, and refractive indices shown in 
Table 1 (although we always assumed the same refractive index for both the foie and 
coarse modes). The parameters in Table 1 are based upon the range of AERONET clima- 
tologies of Dubovik et al. [2002], Our results for this plethora of distributions are shown 
as the small points Figure 6. The circles correspond to the bimodal parameters used to 
calculate the grid of Figure 5 (boldprint values in Table 1), with the size of the circles 
corresponding to the hne mode median radii. The color of the points and the circles 
indicate the volume fraction of aerosols in the bimodal distribution with radii less than 
0.6 /irn, Vfi ne /V tot ai ■ The black lines approximate isolines of constant C f ine / C tota i. 

We first comment on the color scheme shown in Figure 6, which is coded to represent 
the volume fraction of aerosols with radii less than 0.6 /irn. Note that Vfi ne /V to tai is 
not necessarily equivalent to C f ine / C tota i] Vfi ne /V to t a i is the volume fraction of aerosols in 
the size distribution with r < 0.6 /mi, and C fi ne /C tota i is defined by Equation (3). One 
difference between these two definitions is that the hne and coarse mode aerosol sizes 
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described by Equation (3) can overlap, whereas the volume fractions Vf ine and V to t a i can 
not overlap in size because of the hard cut-off at r = 0.6 /jrn. We chose this nomenclature 
for the color code because real atmospheric aerosols are not necessarily constrained by 
Equation (3), and r = 0.6 /jrn is the fine/coarse mode separation radius recommended by 
AERONET (http:/ / aeronet.gsfc.nasa.gov). Also note that the color transitions (i.e., from 
blue to green to orange, etc.) in Figure 6 are roughly parallel to the black lines. This 
is because all lines of constant Vf ine /V tota i or Cf ine /C to tai calculated from the parameters 
in Table 1 have the same general shape. Hence, conclusions about the size distributions 
represented by the circles in Figure 6 can be generalized to all size distributions in Table 
1 . 

Next, we point the reader to the monomodal size distributions in Figure 6 (red circles 
to the right, with C f ine /C tot ai = 1)- In this case, we see that the curvature is always 
negative (02 < 0) and that it decreases as particle size increases; we found this to be true 
for all cases with Cfi ne /C to tai = 1 in Table 1. Positive curvature requires the presence of 
enough coarse mode aerosols to reduce the wavelength dependence of aerosol extinction 
at the longer wavelengths (recall Figure 3). Note that negative curvature can be achieved 
even with the presence of a significant coarse mode component, as evidence by the blue 
points below the 02 = 0 line. 

Also note that the absolute value of the coefficient a\ decreases with increasing particle 
size for the hne monomodal aerosols with C fi ne / C tota i = 1, much like its near cousin a 
(recall that the rightmost red circles in Figure 6 are calculated with the same nine size 
distributions as the uppermost line in Figure 5). This is not the case when significant 
coarse mode aerosols are present and the curvature is positive, as shown by the Cfi ne /C to tai 


DRAFT 


October 11, 2005, 9:45am 


DRAFT 



X- 16 


SCHUSTER ET AL.: ANGSTROM EXPONENT, BIMODAL AEROSOL SIZES 


isolines when a -2 > 0 in Figure 6. This is because positive curvature indicates that the 
coarse mode is already significantly reducing the Angstrom exponent, and increasing the 
effective radius of the fine mode under these conditions increases the radiative efficiency 
of that mode. Hence, the coefficient ci\ shifts from the low magnitudes that are associated 
with coarse mode aerosols to the large magnitudes that are associated with fine mode 
aerosols along the Cfi ne /C to tai isolines with 02 > 0. Finally, we reiterate that the curvature 
exhibited in Figure 6 is always negative for the size distributions of Table 1 when the coarse 
mode is not present. 

The situation where 02 = 0 corresponds to a special case without curvature, and a± = 
—a. It is sometimes postulated that aerosol size distributions without curvature are 
“Junge” distributions (Equation 2), but the large number of points on or near the zero 
line of <22 in Figure 6 indicate that it is possible to have bimodal lognormal aerosol size 
distributions without curvature. The size of the circles and color of the points on the 
line (X 2 = 0 in Figure 6 are also consistent with conventional wisdom about the Angstrom 
exponent (i.e., they indicate that the fine mode median radius decreases and the fine mode 
volume fraction increases as a j increases). 

Note that «i can be significantly different from a, and that the information content of 
sunphotometry measurements can be enhanced by considering the spectral curvature of 
extinction. For instance, recall the demonstration of Figure 3, where extinction calcula- 
tions for two plausible aerosol size distributions produced the same Angstrom exponent. 
One size distribution was monomodal with a fine mode median radius of Rfi ne = 0.21 and 
width of <7 = 0.38, and the other distribution was bimodal with a fine mode fraction of 
C fine / C total = 0.6, a f ine = 0.38, and a crs = 0.75. Both of these size distributions produce 
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an Angstrom exponent of a = 2.0, but they produce measurably different coefficients 
when curvature is considered. Both size distributions can be located by the appropriate 
black squares in Figure 6, revealing that (cq, a 2 ) = (—2.78, —0.76) for the monomodal 
distribution, and (oi, a 2 ) = (—1.74, +0.26) for the bimodal distribution. 

To a close approximation, a = a 2 — oq; hence, we have included two dashed lines in Figure 
6 that represent the traditional guidelines for the Angstrom exponent. Lines of constant 
a run parallel to these lines (not shown). Figure 6 indicates that a 2 — cq > 2 corresponds 
to size distributions dominated by fine mode aerosols and a 2 — cq < 1 corresponds to 
size distributions dominated by coarse mode aerosols, as expected. Intermediate values 
of o 2 — a i (or a) correspond to a wide range of fine mode fractions. 

6. A word about absorption 

The computations shown in Figure 6 represent a single imaginary refractive index of 
m.i = 0.003, or a single bulk absorption coefficient. Natural aerosols have a variety of 
bulk absorption coefficients, however, so we demonstrate the sensitivity of the Angstrom 
exponent to a range of imaginary refractive indices in Figure 7. We used rn t = 10 -6 
to simulate aerosols with negligible absorption and m, = 0.02 to simulate extremely 
absorbing aerosols (in addition to our previous value of rn l = 0.003). We show results for 
three lognormal foie mode fractions of = 0.1, 0.3, and 0.8. 

Figure 7 indicates that the Angstrom exponent is not sensitive to the bulk absorption 
coefficient for the two lowest lognormal foie volume fractions (i.e., c / ine = 0.1 and 0.3). 
We can understand this by noting that a < 1 at these volume fractions, and recalling that 
the spectral absorption of small particles is often approximated as r abs oc A ^ 1 [ Bergstrom , 
2002], Since the spectral dependence of absorption is similar to the spectral dependence 
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of extinction for coarse mode particles, the bulk absorption coefficient has little effect on 
the extinction exponent of these particles. However, the bulk absorption coefficient does 
have an impact on the Angstrom exponent for large volume fractions of fine mode aerosols 


(re., 


Cf ine 

Ctotal 


= 0.8 in Figure 7). Note that the sensitivity is small, though, as the extreme 


variability of rn t in Figure 6 produces a maximum change in a of 0.3 at ^ ine ) = 0.8. 


7. Implications for measurements 

Thus far, we have only discussed the calculated spectral dependence of aerosol extinction 
for hypothetical size distributions of idealistic spherical aerosols. The question arises: 
How does this information convey into real world measurements? Fortunately, we have 
the Aerosol Robotics Network (AERONET) of narrow held of view radiometers, which 
provides aerosol optical depth measurements at over 180 locations worldwide [ Holben 
et al, 1998, 2001]. Additionally, the AERONET radiometers periodically perform sky 
radiance scans which enable the retrieval of column-averaged aerosol size distributions 
and refractive indices [ Dubovik and King, 2000; Dubovik et al, 2000]. In this section, 
we use data from the Aerosol Robotics Network (AERONET) of surface radiometers to 
further illustrate the relationships between the Angstrom exponent, curvature, fine mode 
volume fraction, and fine mode effective radius. 

We begin our discussion with Figure 8, which shows the fine mode volume fractions 
(upper panel) and fine mode effective radii (lower panel) obtained from AERONET size 
distribution retrievals at 53 locations in the years 2000 and 2001 as a function of the 
Angstrom exponent. The AERONET size distributions are provided as column concen- 
trations at 22 radii and not necessarily constrained by Equation (3), so we use the nomen- 
clature Vfine/Vtatai to denote the volume fraction of aerosols with radii less than 0.6 /mi. 
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We calculated the Angstrom exponents and curvatures in Figure 8 from multi-wavelength 
AOT measurements (A = 0.34,0.38,0.44,0.5,0.67,0.87, and 1.02 jum) using least squares 
fits of Equations (6) and (7). We only considered AERONET retrievals with aerosol op- 
tical thicknesses greater than 0.2 at the 0.44 /jm wavelength and a real refractive index 
difference of less than 0.1 between the 0.44 and 1.020 /jrn wavelengths. We also present 
the Angstrom exponent calculations used in Figure 5 as solid lines in the lower panel of 
Figure 8 (note that the variables in the lower panel of Figure 8 are the same as Figure 5, 
with the ordinate and abscissa interchanged). 

Momentarily ignoring the color scheme in Figure 8, we see that the multi- wavelength 
Angstrom exponent provides at least some information about the volume fraction of 
aerosols in the fine mode, but little or no information about the fine mode effective 
radius. For instance, the upper panel of Figure 8 indicates Angstrom exponents greater 
than 2 correspond to fine mode volume fractions greater than about 0.5, Angstrom expo- 
nents less than 1 correspond to fine mode volume fractions less than about 0.5, and that 
intermediate values for a correspond to fine volume fractions of 0.2-0.85. This is consis- 
tent with conventional wisdom that size distributions with large Angstrom exponents are 
dominated by fine mode particle sizes. However, the lower panel of Figure 8 indicates no 
relationship between the Angstrom exponent and the fine mode effective radius. 

The color scheme of Figure 8 indicates that including spectral curvature in our analysis 
may enhance our knowledge about the volume fraction and effective radius of fine mode 
aerosols at intermediate values of the Angstrom exponent (1 ot < 2), which we discuss 
here. Suppose a value of a = 1.5 was calculated from multi-wavelength AOTs, which we 
have already noted corresponds to a fine mode volume fraction range of 0.2-0.85. Paying 
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attention to the color scheme in Figure 8, we see that the negative curvatures favor high 
volume fractions and positive curvatures favor low volume fractions at this intermediate 
value for a. Our analysis of Figure 8 at a — 1.5 indicates that a 2 < —0.3 (blue squares) 
corresponds to fine mode volume fractions greater than 0.4. Similarly, a 2 > 0.3 (red 
squares) corresponds to fine mode volume fractions less than 0.65. Likewise, the largest 
fine mode effective radii require large negative curvatures (lower panel of Figure 8), con- 
sistent with our theoretical discussion of Section 5. Hence, the presence of curvature can 
be used to improve our assessment of aerosol size distributions at intermediate Angstrom 
exponents. Note that the curvature adds no new information for a < 1 and a > 2. 

We plot the same AERONET data in Figure 9, this time mimicking the format of 
Figure 6. The data indicate that the coefficient cq and the curvature a 2 tend to decrease 
as the fraction of particles in the fine mode increases, consistent with our theoretical plot 
of Figure 6. This is consistent with Eck et al. [1999]; Eck. et al. [2001b], and Kaufman 
[1993], who found that the foie mode dominates when the curvature is negative, and that 
the coarse mode contributes significantly when positive curvature is present. Nonetheless, 
the data also indicate that negative curvatures are possible for size distributions dominated 
by coarse mode particles (i.e., blue points below the a 2 = 0 line) and positive curvatures 
are possible for size distributions dominated by fine mode particles (red points above the 
a 2 = 0 line). 

In satellite remote sensing, the retrieval of the Angstrom exponent is based upon a 
single pair of wavelengths [ Nakajima and Higurashi, 1998; Higurashi and Nakajima, 1999; 
Deuze et al., 2000; Ignatov and Stowe, 2002; Jeong et al., 2005]. Conclusions about 
the size distribution in such cases must consider the retrieval wavelengths, as different 
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wavelength pairs will provide different information about the aerosol size distribution. 
This is demonstrated in Figure 10, where we have plotted the Angstrom exponents of 
two-wavelength pairs for synthetic bimodal lognormal distributions and for AERONET 
data. Here, the AERONET data is once again obtained from almucantar retrievals at 
53 AERONET sites, and the circles are calculated from Equations (3) and (4) using the 
boldprint bimodal parameters of Table 1. The Angstrom exponents for each axis were 
calculated using only two wavelengths in Equation (6). 

Note that the Angstrom exponent from the AERONET dataset is rarely linear, as 
evidenced by the small fraction of data on the 1:1 line in Figure 10. The upper panel 
indicates that negative curvature (<22 < 0) favors larger fractions of fine mode aerosols, but 
fine mode volume fractions greater than 0.5 are not required for negative curvature. Like- 
wise, the lower panel of Figure 10 indicates that positive curvature (<22 > 0) favors smaller 
fine mode radii, but this is not always the case. Separation of colors for AERONET data 
would indicate that one could obtain fine volume fraction and effective radius from this 
type of plot, but there is considerable intermingling of colors, making detailed conclusions 
about the size distribution nebulous. 

This is illustrated more clearly in Figure 11, where we present the same AERONET 
data once again. Figure 11a indicates that calculating the Angstrom exponent based 
upon the 0.38 and 0.44 /mi wavelengths will tell us very little about the fine mode volume 
fraction of aerosols, as evidenced by the large amount of scatter and the small slope of the 
linear regression. This is because both of these wavelengths are less sensitive to coarse 
mode aerosols than the 0.5 /irn wavelength shown in Figure 1. This lack of sensitivity to 
the coarse mode increases the sensitivity of these wavelengths to changes in the fine mode 
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effective radius, as shown in Figure lib. Conversely, the Angstrom exponents based upon 
the 0.67 and 0.87 /jrn wavelengths are indeed sensitive to the fine mode volume fraction 
(albeit with much scatter), as shown by the slope in Figure 11c, and have nearly no 
sensitivity to the fine mode effective radius (Figure lid). This wavelength sensitivity of a 
is consistent with the work of others [Eck. et al, 2001b; O'Neill et al., 2001a; Reid et al., 
1999], 

8. Real versus synthetic data 

Note that the color transitions for the AERONET data in Figure 9 are not nearly as 
clean as the transitions for the synthetic data of Figure 6, which we explain in this section. 
Measured AOTs were used to calculate the coefficients (ai,a 2 ) in Figure 9, so this is our 
benchmark. Nonetheless, measurement noise exists in Figure 9 that is not considered in 
the synthetic calculations of Figure 6. Sunphotometry measurements of AOT are typically 
accurate to 0.02 or better for newly calibrated field instruments [ Holben et al., 1998], and 
this uncertainty alters the Angstrom exponent a by 0.03-0.04; we can expect similar 
variability in a\. This variability in a\ (and a 2 ) manifests itself as an intermingling of 
colors in Figure 9, as does any errors associated with the AERONET size distributions 
retrievals. The data in Figure 6, on the other hand, are somewhat idealistic because these 
noise issues are not taken into consideration. We chose a single imaginary refractive index 
for the computations of Figure 6, which reduces the computational noise associated with 
the discussion of Section 5 

We also chose to use y 2 iteration to infer aerosol absorption for Figure 6, which gives 
a spectrally smooth and plausible refractive index for an internal mixture of ammonium 
sulfate, water, and black carbon. We note that using y 2 iteration of the refractive index 
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with AERONET distributions (which are not necessarily lognormal) also produces an 
idealistic scenario similar to Figure 6 (not shown), so the shape of the size distributions 
is not a source of noise. We also note that we can reproduce measured AOTs to within 
0.03 using the AERONET size distributions and AERONET refractive indices at the four 
almucantar wavelengths, so our Mie computations are correct. 

Other issues besides the measurement noise are also responsible for the differences be- 
tween Figures 6 and 9. Applying the y 2 iteration procedure to the AERONET aerosol re- 
trievals does not necessarily produce refractive indices that exactly match the AERONET 
refractive indices (although the magnitude of the difference is limited by the magnitude 
of the spectral variability of the AERONET retrieval). Aerosol optical depth calculations 
are very sensitive to the real refractive index, as variations in the real refractive index of 
0.02 at the 0.440 /ini wavelength can produce optical depth variations as large as 0.12. 
Spectral variations of the refractive index beyond the small spectral variability for water 
and ammonium sulfate will not show up in Figure 6, so some of the noise in Figure 9 could 
occur because of realistic spectral variations in the AERONET refractive index. Large 
spectral variability of the real refractive index is not expected for nondust aerosols, how- 
ever. Cases where the retrieved real refractive index exhibits unrealistically large spectral 
variability may also indicate an erroneous fine mode fraction, since the AERONET re- 
trieval is constrained by the AOT at the four almucantar wavelengths. Such cases would 
also appear as noise in Figure 9. 

9. Conclusions 

We have discussed the relationship between polydisperse aerosol size distributions and 
the spectral dependence of the aerosol optical thickness, using both Mie calculations and 
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AERONET aerosol retrievals. We began by showing the expected inverse relationship 
between the Angstrom exponent and the effective radius for synthetic monomodal aerosol 
size distributions. 

Next, we discussed the sensitivity of the Angstrom exponent and curvature to the details 
associated with bimodal aerosol size distributions, focusing specifically on the fine mode 
effective radius and the fine mode fraction of aerosols. We demonstrated that calculations 
of the Angstrom exponent for seven wavelengths (0.34, 0.38, 0.44, 0.50, 0.67, 0.87, and 
1.02 /ini) are sensitive to the fraction of aerosols in the fine mode ( Vf ine /V tota i ), but 
not to the fine mode effective radius (. R e ff(fine )). Nonetheless, this multi- wavelength 
Angstrom exponent is not a rigorous indicator of Vfi ne . Rather, AERONET retrievals at 
approximately 50 sites indicate that a > 2 correspond to fine mode fractions of Vf ine > 0.5, 
and a < 1 correspond to fine mode fractions of Vj ine < 0.5. The same AERONET dataset 
indicates that intermediate values of a correspond to fine mode fractions between 0.2 
and 0.85, but improved estimates of the fine mode aerosol fractions can be obtained by 
considering the spectral curvature. For instance, if a = 1.5, then curvatures of a 2 < —0.3 
indicate fine mode fractions of 0.5 < Vf ine < 0.85, and a 2 > +0.3 indicate fine mode 
fractions of 0.2 < Vf ine < 0.5. Although the curvature is useful for estimating the fine 
mode fraction of aerosols at intermediate values of a, it provides no additional information 
when a < 1 and a > 2. 

Spectral curvature of aerosol extinction also plays a significant role when the Angstrom 
exponent is calculated with only two wavelengths. Angstrom exponents calculated from 
longer wavelength pairs (A = 670, 870 //m) are sensitive to the fine mode fraction of 
aerosols but not the fine mode effective radius; conversely, shorter wavelength pairs (A = 
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380,440 /jrn) are sensitive to the fine mode effective radius but not the fine mode fraction 
(see also Eck et al. [1999]; Eck. et al. [2001b]; Kaufm.an [1993]). Hence, it is important 
to consider the wavelength pair used to calculate the Angstrom exponent when making 
qualitative assessments about the corresponding aerosol size distributions. 
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Figure 1 . Relative optical thickness for monomodal distributions of aerosols described 
by Equation (3) with constant C i, C 2 = 0, ay = 0.38, and variable effective radii. 
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Figure 2. Relationship between the effective radius and the modal median radius of 
Equation (3). The discrepancy between effective radius and modal median radius increases 
as the mode width increases. 
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Figure 3. Spectral extinction for two different size distributions with the same Angstrom 
exponent (a = 2). The squares correspond to a monomodal lognormal size distribution 
with Rfi ne = 0.21 and a = 0.38; the circles correspond to a bimodal lognormal size 
distribution with Cf ine /C t otai = 0.6, R fine = 0.12 , R crs = 3.2 ,a fine = 0.38, and a crs = 
0.75. A refractive index of m = 1.37 — 0.003f was used in both cases. 
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Figure 4. The Angstrom exponent (upper panel) and curvature (lower panel) 
for monomodal lognormal aerosol size distributions with varying effective radii. The 
Angstrom exponent is sensitive to fine monomodal aerosol size distributions {R e ff < 
0.25 /mi), but not to coarse monomodal aerosol size distributions (R e ff > 2/un). The 
curvature is negative for fine monomodal aerosols. 
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Fine mode effective radius (pm) 

Figure 5. Calculated Angstrom exponents for bimodal aerosol size distributions using 
five fine mode fractions and nine fine mode median radii in Equation (3). The coarse 
mode radius was held constant at R crs = 3.2 //m and the mode widths were held constant 
at a fi ne = 0.38, a crs = 0.75, and the refractive index is m = 1.37 — 0.003?' for both modes. 
These parameters are listed in boldprint in Table 1. 
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Figure 6. Volume fraction of fine mode aerosols as a function of the coefficients a± 
and 02 in Equation (7) for 13,860 bimodal aerosol size distributions parameterized with 
the values in Table 1. The circles correspond to the boldprint values in Table 1 and 
are sized relative to the fine mode radius. The black squares correspond to the two size 
distributions with a = 2 used in Figure 3. The dashed lines closely approximate lines of 
constant a. 
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Figure 7. Sensitivity of Angstrom exponent to the aerosol imaginary refractive index 

(or equivalently, the bulk absorption coefficient) for bimodal aerosol distributions with 

three different fine volume fractions. The Angstrom exponent shows some sensitivity 

to absorption when the bimodal size distribution is dominated by fine mode aerosols 
c 

( c f t ™ e i = 0.8), but shows virtually no sensitivity to absorption at the smaller fine mode 
aerosol fractions. 
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Figure 8. Fine volume fraction and fine mode effective radius as a function of the 
multi-wavelength Angstrom exponent. Squares correspond to AERONET almucantar 
retrievals and coincident AOTs at 53 locations in the years 2000 and 2001; the color of 
the squares correspond to the spectral curvature of the sunphotometry measurements, 
and the size of the squares correspond to the magnitude of the curvature, |a 2 |. Lines in 
the lower panel represent constant C f ine / C tota i values for bimodal lognormal aerosol size 
distributions parameterized with the boldprint values in Table 1. 
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Figure 9. Same as Figure 6, except that the small points indicate the fine volume 
fractions obtained from AERONET retrievals. 
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Figure 10. Fine mode volume fraction (upper panel) and fine mode effective radius 
(lower panel) as a function of two- wavelength Angstrom exponents. Small points are 
sunphotometry data from 53 AERONET sites. Circles are synthetic lognormals with 
boldprint values in Table 1, with sizes proportional to the fine mode median radius. 
Shorter wavelength Angstrom exponents are more sensitive to fine mode particle size 
than longer wavelength Angstrom exponents; longer wavelength Angstrom exponents are 
more sensitive to fine volume fraction than shorter wavelength Angstrom exponents. 
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Figure 11. Fine mode volume fraction and fine mode effective radius as a function of the 
380-440 two- wavelength Angstrom exponent (left panels) and the 670-870 two-wavelength 
Angstrom exponent (right panels). A comparison of the linear regression slopes in each 
panel indicates that the long wavelength Angstrom exponent has greater sensitivity to the 
fine mode aerosol fraction than the short wavelength Angstrom exponent (upper panels). 
Likewise, the short wavelength Angstrom exponent has greater sensitivity to the fine mode 
effective radius than the long wavelength Angstrom exponent (lower panels)). 
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Table 1. Bimodal lognormal size distribution parameters used to calculate the points in 
Figures 6 and 10. Values in boldprint are also used to calculate the circles in Figures 5-10. 


Parameter 

Value 

^/me(^ 0 

0.06, 0.09, 0.12, 0.15, 0.18, 0.21, 0.24, 0.27, 0.30 

O fine 

0.38, 0.5 

F^crs 

1.9, 2.2, 2.7, 3.2, 3.7 

&crs 

0.75, 1.0 

m r 

1.34, 1.37, 1.4, 1.43, 1.47, 1.5, 1.54 

rrii 

0.003 

C fine/ Cfotal 

0.0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0 
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